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ABSTRACT 

We investigate chemical abundance variations along the sightlines through 13 damped Lya systems 
(DLA). We introduce a technique designed to identify abundance variations in multiple velocity bins 
along the sightlines and perform a series of Monte-Carlo simulations to derive quantitative limits and 
values. The majority of these DLA have very uniform relative abundances: 95% c.l. upper limits of 
0.2 dex dispersion along the sightline and best fit values typically < 0.1 dex. This surprising result 
indicates the gas comprising individual DLA has similar enrichment history, nearly identical differential 
depletion, and small abundance variations relating to ionization corrections. This uniformity contrasts 
with stellar abundance variations observed within the Calaxy and with the differences in abundances 
observed between galaxies of the Local group. It also contrasts with variations in differential depletion 
along sightlines through the SMC, LMC, and Milky Way. The results constrain the process of metal 
production and dust formation in the early universe and reflect on processes of galaxy formation. In 
terms of depletion the observations indicate a very low porosity of significantly depleted gas, substantially 
lower than the filling factor observed in present-day metal-poor galaxies. Finally, the observed chemical 
uniformity may present a difficult challenge to scenarios which assume individual DLA are comprised of 
multiple protogalactic clumps. 

Subject headings: galaxies: formation, galaxies: high-redshift, quasars: absorption lines 



1. INTRODUCTION 

Studies of chemical abundance ratios in Calactic metal- 
poor stars address processes of nucleosynthesis, star forma- 
tion, and ultimately galaxy formation in the early universe. 
Of particular importance is revealing variations or trends 
in the relative abundances which identify distinct stellar 
populations, unique nucleosynthetic sites, and/or provide 
insight into the time-scales and channels of star formation. 
In very low metallicity ([Fe/H] < —3) Galactic stars, for 
example, observers have identified important trends in the 
Fe-peak elements and large variations in several abundance 
ratios (e.g. Co/Fe, Sr/Fe) which suggest unique nucleosyn- 
thetic processes and inhomogencous chemical enrichment 
(e.g. McWiUiam et al. 1995; Johnson 2002). In contrast, 
the abundance patterns of the majority of Galactic stars 
with —2 < [Fe/H] < —1 are remarkably similar (Gratton 
& Sneden 1991; Fulbright 2000), suggesting a rapid en- 
richment phase for the Galactic halo over this metallicity 
range. At higher metallicity, there are important trends in 
the relative abundances, in particular the steady decline 
of the a-elements (0,S,Si) to the Fe-peak, which mark the 
onset of Type la SN enrichment (Tinsley 1979). Compar- 
isons of these relative abundances with extragalactic stel- 
lar populations help map paths of galaxy formation (e.g. 
Matteucci & Recchi 2001). In the few galaxies where 
detailed studies can be performed (the SMC and dwarf 
spheroidals; Venn 1998; Shetrone, Cote, & Sargent 2001; 
Bonifacio et al. 2000), one often observes significant differ- 
ences with stars at the same metallicity in the Milky Way. 
Future studies of entire stellar populations within these 
galaxies will help refine their star formation histories. 

^ Visiting Astronomer, W.M. Keck Telescope. The Keck Observa- 
tory is a joint facility of the University of California and the Califor- 
nia Institute of Technology. 



An examination of relative abundances of gas in lo- 
cal galaxies complements these stellar abundance stud- 
ies. Sightlines through the Galactic ISM probe gas with 
a range of physical conditions (e.g. volume density, dust- 
to-gas ratio, ionization state). Therefore, the observed 
gas-phase abundance ratios along a sightline may vary by 
an order of magnitude as the sightline penetrates clouds 
arising in various phases of the ISM (e.g. Savage and Sem- 
bach 1996). Differences in nucleosynthetic enrichment or 
photoionization are likely to be small such that the varia- 
tions primarily reflect differences in depletion level or dust 
composition. The LMC and SMC also exhibit variations 
in these gas-phase abundance ratios (Welty et al. 1999, 
2001) and together the observations impact models for the 
formation and distribution of dust as well as the physical 
state of galactic interstellar media. 

Observations of the damped Lya systems (DLA) - quasar 
absorption line systems with A^(HI) > 2 x 10^*^ cm^^ - 
probe the ISM of high z galaxies. Surveys of the DLA 
reveal chemical evolution in the early universe (Pettini et 
al. 1994, 1997; Prochaska & Wolfe 2000), and exam- 
ine the dust depletion, nucleosynthetic enrichment, and 
ionization of these protogalaxies (Lu et al. 1996; Vladilo 
1998; Prochaska & Wolfe 2002; Prochaska et al. 2002a). 
Prochaska & Wolfe (2002) recently reviewed the chem- 
ical abundances of the DLA and emphasized the unifor- 
mity in the relative abundance patterns (e.g. Si/Fe, Ni/Fe, 
Zn/Si) of their « 30 systems. Even though damped Lya 
metallicities span over an order of magnitude, the systems 
show very similar abundance ratios suggesting protogalax- 
ies have common enrichment histories. The authors do 
identify a mild trend of increasing Si/Fe and Zn/Fe with 
increasing metallicity which is well explained by the ef- 
fects of dust depletion, yet even this trend spans only a 
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factor of w 3 in the ratios. In addition to the uniformity 
in abundance patterns from galaxy to galaxy, Prochaska 
& Wolfe (2002) asserted that this uniformity holds along 
the sightlines penetrating each protogalaxy. That is, the 
individual 'clouds' comprising mctal-linc profiles show the 
same relative abundances within a given damped Lya sys- 
tem. If verified, this would have important implications 
for the ISM of high z galaxies and the enrichment of gas 
in the early universe. 

Prochaska & Wolfe (1996) first quantitatively investi- 
gated variations in the chemical abundances of a single 
damped Lya system. They compared Zn, Si, Fe, Ni, and 
Cr ionic column densities along the observed velocity pro- 
files using the apparent optical depth method (Savage and 
Sembach 1991). Their analysis showed that the chemi- 
cal abundances were uniform to within statistical uncer- 
tainties. They then argued that the absence of significant 
variations demonstrates a low depletion level within this 
protogalaxy. More recently, Lopez et al. (2002) performed 
the first detailed line-profile analysis of the relative abun- 
dances for a different z > 2 damped system and also found 
nearly constant column density ratios among the compo- 
nents comprising the line-profile solution (see also Petit- 
jean, Srianand, & Ledoux 2002). Together, these studies 
argue the gas within high z protogalaxies has a similar 
enrichment history and uniform differential depletion. 

In this paper, we search for variations in the relative 
abundances along the sightlines of 13 damped Lya sys- 
tems. These systems were selected to have high signal- 
to-noise velocity profiles with velocity widths exceeding 
40 kms""'^. The latter criterion ensures a large enough 
baseline to examine variations within the protogalaxy. We 
introduce a quantitative method which can be applied to 
our observations as well as mock spectra derived from nu- 
merical simulations and semi-analytic models of galaxy 
formation. The technique assesses changes in the phys- 
ical conditions and enrichment histories of the gas within 
each damped system. In turn, our results place valuable 
constraints on the various morphological models of the 
damped Lya systems and, ultimately, scenarios of galaxy 
formation. 

2. ANALYSIS OF THE OBSERVATIONS 

To investigate variations in the relative abundances of 
the damped systems, we focus on pairs of low-ion transi- 
tions for 13 damped systems from our damped Lya database^ 
(Prochaska ct al. 2001). We selected damped Lya systems 
with transition pairs that satisfied the following criteria; 
(1) signal-to-noise (S/N) > 15 per pixel; (2) a velocity 
width Av > 40 kms~^ to provide a base-line for exam- 
ining abundance variations; (3) transitions which are not 
saturated but have peak optical depth with > 5(t signifi- 
cance; (4) a transition pair with different refractory char- 
acteristics, nucleosynthetic origin, and/or photoionization 
dependencies. With respect to the last criterion, the a- 
ions (e.g. Si"*", S"*") are ideal complements to the Fe-peak 
ions (e.g. Ni+, Fc+). The former are non-refractory or 
mildly refractory, are primarily produced in Type II SN 
and are less ionized in HII regions, while the latter are 
refractory, primarily produced in Type la SN, and in the 
case of Fe"*" are more highly ionized. Table 1 summarizes 

^ http://kinpgin.ucsd.edu/~hiresdla/ 



the damped Lya systems examined in this paper and lists 
the transitions comprising our analysis. In nearly every 
case, we observe multiple transitions for a given damped 
Lya system and are confident that the transitions listed 
in Table 1 are free of line-blending. 

We considered two techniques for probing variations in 
the relative abundances of a pair of metal-line profiles: (1) 
a comparison of the ionic column densities derived from 
the apparent optical depth (Savage and Sembach 1991) 
in a series of velocity bins: and (2) a comparison of the 
column densities derived from Voigt profile fitting of in- 
dividual velocity components. The latter approach has 
the advantage that one presumably isolates abundances 
in specific physical components (e.g. 'clouds' or 'clumps'), 
whereas the velocity bins may contain one or more such 
components. There are several weaknesses, however, to 
the Voigt profile analysis. First, identifying the component 
structure of a metal-line profile is both time consuming 
and somewhat arbitrary. One typically introduces enough 
'clouds' to produce a line-profile solution that yields a re- 
duced xt ~ 1- Therefore, high S/N profiles invariably 
exhibit more components than lower quality data. Sec- 
ond, blended components generally exhibit degeneracies 
between their column densities. In order to accurately 
compare ionic column densities one must bin these blended 
components and carefully consider their correlated errors. 
Third and most important, in cases where a pair of tran- 
sitions have very large abundance differences it may be 
impossible to derive a self-consistent profile solution. For 
example, one ion may show several components which are 
not observed in the other ion. In these cases, one cannot 
compare the abundances component by component. The 
risk is that this would bias the results against systems with 
large abundance variations. 

In this paper, we concentrate on the apparent optical 
depth method. For many of the sightlines presented here, 
we have also examined the results from a Voigt profile 
analysis and find similar results (e.g. the systems included 
in an analysis of the fine structure constant; Murphy et al. 
2001). After experimenting with the size of velocity bins, 
we found that binning the apparent optical depth over 
10 pixels (20 kms~^) offered the best comprise between 
minimizing pixel-to-pixel noise and maximizing velocity 
resolution. For the very highest S/N data, we reduced the 
10 pixel bin to 5 pixels. In all cases, we clipped velocity 
bins with statistical significance less than 2.5cr, i.e., in re- 
gions where there is very little absorption. At 20 kms~^ or 
even 10 kms~^, it is possible that the velocity bins include 
two or more 'clouds' with differing abundances which are 
averaged to a single value; we do not resolve abundance 
variations on this small scale. On the other hand, the 
typical separation between the main components compris- 
ing the metal-line profiles is > 20kms~^. Our analysis 
focuses on these main components. 

Figure 1 presents the results for the damped Lya system 
at 2 = 1.92 toward Q2206-19 which is roughly represen- 
tative of the full sample (Figure 2). The upper panels 
show the Sill 1808 and Fell 1611 velocity profiles while 
the lower panel plots the ratio logj^Q (Si^/Fe"*") normalized 
to solar relative abundances for the statistically significant 
velocity bins. The dashed-line in the lower panel indicates 
the ratio of the total column densities derived from these 
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Table 1 
OBSERVATIONAL SAMPLE 



QSO 






A?; 


Transition 


[X/Y]t 






rr b 


xi 




^all 


^best 








(kms-1) 


Pairs 


(dex) 


(kms-1) 




(dex) 




(dex) 


(dex) 


(dex) 


PH957 


2.309 


-1.46 


54 


Znll 2026 
Nill 1741 


0.22 


10 


3 


0.18 


1.65 


0.24 


0.47 


0.03 


Q0201+36 


2.463 


-0.41 


202 


Sill 1808 
Fell 1608 


0.45 


20 


12 


0.07 


3.21 


0.38 


0.15 


0.09 


Q0347-38 


3.025 


-1.17 


94 


Sill 1808 
Fell 1608 


0.50 


10 


5 


0.47 


11.0 


0.74 


0.90 


0.39 


Q0458-02 


2.040 


-1.19 


86 


Znll 2026 
Cr II 2056 


0.39 


20 


4 


0.48 


4.81 


0.35 


0.45 


0.15 


HS0741+47 


3.017 


-1.69 


40 


S II 1259 
Fell 1608 


0.24 


10 


5 


0.08 


2.23 


0.23 


0.22 


0.07 


PSS0957+33 


4.180 


-1.50 


346 


Sill 1304 
Fell 1608 


0.37 


20 


6 


0.26 


5.34 


0.70 


0.68 


0.29 


Q1223+17 


2.466 


-1.59 


94 


Sill 1808 
Nill 1751 


0.13 


10 


10 


0.11 


1.82 


0.40 


0.19 


0.09 


Q1331-17 


1.776 


-1.30 


70 


Zn II 2026 
Fell 2374 


0.82 


10 


9 


0.23 


22.8 


0.66 


0.41 


0.24 


GB1759+75 


2.625 


-0.82 


78 


Sill 1808 
Nill 1751 


0.35 


10 


8 


0.18 


2.50 


0.44 


0.27 


0.11 


Q2206-19 


1.920 


-0.42 


130 


Sill 1808 
Fell 1611 


0.40 


20 


7 


0.06 


1.17 


0.24 


0.15 


0.04 


Q2230+02 


1.864 


-0.74 


177 


Sill 1808 
Nill 1741 


0.23 


20 


7 


0.08 


1.30 


0.31 


0.18 


0.03 


Q2231-00 


2.066 


-0.86 


127 


Sill 1808 
Nill 1741 


0.28 


20 


4 


0.08 


0.93 


0.26 


0.33 


0.03 


Q2359-02 


2.154 


-1.58 


127 


Sill 1526 
Fell 1608 


0.32 


20 


4 


0.08 


0.97 


0.24 


0.32 


0.04 



"Metallicity derived from Zn, Si, or S 
''Logarithmic RMS dispersion in X/Y 

Note. — Asngh ^alh ^best are measures of the intrinsic abundance variation along ea^h sightline. See § 3 
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Fig. 1. — The top panels show the velocity profiles for the Sill 1808 and Fell 1601 transitions for the damped Lya system at 2 = 1.92 
toward Q2206-19. The lower panel plots the logarithmic ratio of the column densities of the two ions measured using the apparent optical 
depth method in 20 kms~^ bins where the dashed line indicates the total logarithmic ratio. Deviations from this total value are small for the 
7 bins (u < 15%) and they are consistent with expectation from statistical error alone. 



two profiles. We compare the individual bins against this 
total ratio to examine departures from uniformity. The 
error bars only reflect statistical error; we do not account 
for systematic errors like continuum placement. In the 
systems analysed in this paper, which have relatively high 
quality data, we expect such systematic errors are small. 
The damped Lya system presented in Figure 1 shows very 
small variations over the entire profile. The scatter in the 
relative abundances is < 15% (w 0.05 dex) of the total ra- 
tio, consistent with the expected statistical error, i.e., the 
reduced ~ 1- 

We have repeated this analysis for the remaining damped 
systems and summarize the results for all of the velocity 
bins for the 13 DLA in Figure 3. The top panel shows the 
departure from uniformity in numbers of standard devia- 
tion n„ where a is derived with standard propagation of 
statistical error. The lower panel presents the departures 
on a logarithmic scale relative to the total X/Y ratio. To- 
gether the two panels describe the magnitude and signifi- 
cance of the variations in our sample. For both panels, the 



X-axis indicates the total enhancement of X/Y relative to 
solar on a logarithmic scale [X/Y]t- Higher [X/Y]t val- 
ues indicate systems with larger nucleosynthetic enhance- 
ments, dust depletion, and/or photoionization levels. 

Consistent with a qualitative inspection of the damped 
Lya profiles, there are only two systems with deviations 
> 0.3 dex and only a few others have statistically sig- 
nificant deviations at greater than 0.2 dex. Of the two 
systems with the largest departures from uniformity, one 
is the damped system at z = 1.77 toward Q 133 1-1-1 7 which 
shows the largest [X/Y] value and might have been most 
expected to show significant deviations. The other system, 
at z = 3.025 toward Q0347-38, was previously known to 
have abundance differences among its metal-line profiles 
and H2 absorption lines (Levshakov et al. 2002). In- 
specting Figure 3, one infers that the abundance varia- 
tions within most damped systems are at levels less than 
0.2 dex. In the following section, we quantitatively address 
the maximum variations allowed by the entire sample and 
the maximum deviations allowed by individual damped 
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Fig. 2. — Same as Figure 1 but for the remaining damped Lya systems comprising the sample. 
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Fig. 3. — The two panels summarize the (a) statistical significance and (b) logarithmic magnitude of the deviations from uniformity for all 
of the velocity bins of our DLA sample as a function of the integrated X/Y ratio, [X/YJy. Only two systems exhibit variations at greater 
than 0.3 dex and the majority of bins have values less than 0.1 dex. 



systems. 

3. MEASURES OF INTRINSIC VARIATION 

To place quantitative constraints on the nucleosynthetic 
and differential depletion variations in the DLA, we will 
now determine the maximum abundance variations allowed 
by our observations. For a crude measure of the typcial 
abundance variation, we evaluate the maximum component- 
to-component variation allowed when considering the en- 
tire sample of damped Lya systems together, i.e., the 
maximum intrinsic abundance variation which is consis- 
tent with the complete distribution of observed deviations 
(Figure 3). To this end, we place limits on the intrinsic 
abundance dispersion S which is the combined variation 
due to nucleosynthetic enrichment, differential depletion, 
and photoionization effects. Specifically, S defines a uni- 
form distribution of intrinsic abundance variations for the 
Nh velocity bins comprising a single DLA. 

For 6 values ranging from to 0.3 dex, we perform a se- 
ries of Monte-Carlo analyses. In each Monte-Carlo run, we 
simulate 1000 damped systems each with Nb velocity bins. 



For every velocity bin, we randomly draw four quantities: 
(1) the statistical uncertainty a in the abundance ratio, 
drawn from the observed distribution of uncertainties for 
all bins of all of the damped systems in our sample; (2) 
the contribution to the observed ratio from this statisti- 
cal uncertainty assuming a Gaussian distribution A^rror 
with mean of zero and variance cr^; (3) a intrinsic abun- 
dance variation Aatnd drawn uniformly^ from the interval 
Aabnd £ [^<5, 5]] and (4) a weighting factor Wi for each bin 
drawn uniformly from 1 to 10 in order to simulate the 
observed variability in component strengths. We then cal- 
culate the total ratio: 

E ■ 10^' 



^ Note that we assume uniform distributions for the abundance de- 
viations throughout this paper. If we adopted a normal distribution, 
we would place even tighter constraints on the abundance variations. 
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Fk;. 1. Confidence level plot for a scries of Monte-Carlo simulations assuming a range of intrinsic abundance variation <5. The dotted 
lines indicate the approximate 95% c.l. values and the dashed lines refer to 99% c.l. These simulations indicate that the damped Lya systems 
do not allow intrinsic variations greater than 0.25 dex at the 99% c.l. At the same time, the observations do require an intrinsic variation at 
> 0.05 dex. 



where Aj = Aahnd + A^rror which corresponds to the 
[X/Y] — [X/Y]t vahics displayed in Figure 3. For each 
velocity bin in every damped system, we determine the 
statistical significance 

n„ = {T- Ai)/Aerror (2) 

of the variation and for the 1000 trials determine a dis- 
tribution of rirj values. This Monte-Carlo distribution of 
rio- values is then compared against the observed distribu- 
tion (top panel of Figure 3) with a two-sided Kolmogorov- 
Smirnov (K-S) test. The K-S statistical test focuses on 
the median of two distributions and the K-S probabil- 
ity gives the likelihood that two distributions could have 
been drawn from the same parent population (Press et al. 
1992). 

In Figure 4 we present the results of our Monte-Carlo 
simulations as a confidence-level plot for a range of 5 val- 
ues. The 99% c.l. limits for 5 are 0.07 and 0.23 dex. If a 
single abundance process dominates the abundance vari- 
ations, these are the limits on that process. If multiple 
processes contribute (e.g. differential depletion and nucle- 
osynthetic enrichment), then the limits on each process 
are even tighter. Note, these conclusions are insensitive to 
the value of Ni, for Nb > 2. The results in Figure 4 quan- 
titatively support the characteristic of uniformity readily 
observed in the metal-line profiles. If depletion and nucle- 
osynthetic variations arc present in the damped systems, 
they are present at very small levels in the majority of pro- 
togalaxics. While the simulations imply strict upper limits 
to the magnitude of abundance variations, Figure 4 also 
reveals that the observations are inconsistent with zero 



abundance variations. At the 95% c.l., we find that 6 must 
exceed w0.05 dex to match the observations^. Therefore, 
the data do require some intrinsic abundance variations. 

Impressed by the absence of abundance variations in the 
detailed profiles as well as the integrated [X/Y]t values, 
wc posed the following hypothesis: Are the abundance 
ratios within individual damped systems consistent with 
being drawn from a single population described by the 
distribution of [X/Y]t values? To address this hypothe- 
sis, we performed a similar Monte-Carlo simulation where 
the [X/Y] values for each velocity bin arc drawn from the 
observed [X/Y]t for all of the damped Lya systems in 
our database (over 30 DLA; Prochaska et al. 2001), i.e., 
the variations A, are due to differences in [X/Y]t among 
the observed damped Lya systems. Wc simulated 1000 
damped systems with iV(, = 4 bins each and fotmd a KS 
probability Pks = 49% suggesting the dispersion in the 
component-by-component abundances of the DLA is con- 
sistent with the dispersion in the integrated [X/Y]t val- 
ues: The damped systems exhibit a 'scale-free dispersion' 
in their relative abundances. This trait did not have to 
hold for these protogalaxies. One could surmise a scenario 
where the gas 'clouds' comprising a single damped system 
would have the same physical properties and abundances 
ratios but their [X/Y]t values would differ by greater than 
1 dex from system to system. Conversely, the gas compris- 
ing a single damped system could exhibit a range of deple- 

^ It is possible that we are slightly overestimating this lower limit 
(and the upper limits) because we have not included systematic error 
in our analysis, but we expect this would contribute at < 0.05 dex. 
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tion levels and nucleosynthetic enrichment which average 
out to a very similar [X/Y]t ratio in all damped systems. 
In our current sample, however, we find that both the gas 
within a single damped system and the systems together 
have similar dispersion in the relative abundances. 

The results presented in Figure 4 provide an estimate of 
the overall variations in the damped Lya systems. One can 
perform a similar set of Monte-Carlo simulations to inves- 
tigate the deviations allowed within each damped system 
and then examine trends with other physical properties of 
the damped systems (e.g. [X/Y]t, as in Figure 3). We will 
consider three measures of the uniformity within individ- 
ual damped systems. With one measure, wc study an ex- 
treme scenario where all of the variation arises from a sin- 
gle velocity bin. We calculate, Asngi, the minimum [X/Y] 
variation in a single velocity bin which gives Xm-c > Xobs 
in over 95% of the trials in the Monte-Carlo analysis. The 
two other measures assume deviations in all of the veloc- 
ity bins with values drawn from a uniform distribution. 
We calculate (i) Aaii , the minimum variation which when 
applied to every velocity bin gives Xm-c ^ Xobs i'^ over 
95% of the trials; and (ii) Ahest, the variation which ap- 
plied to each component has the highest probability of 
yielding = Xobs ^ 10%- We believe the Aaii values are 
the most realistic upper limits to abundance variations and 
that Abest reflects the most likely value. 

For each damped system, wc take the observed num- 
ber of velocity bins and adopt the measured statistical 
errors. To simplify the analysis, we parameterize all of 
the sources of intrinsic abundance variations into a single 
term as above. Columns 11-13 of Table 1 present the re- 
sults from a series of Monte-Carlo simulations comprised 
of 1000 trials for each damped system. It is important to 
stress that Agngi and Aaii each has a systematic behavior 
related to the number of bins observed. Generally, sys- 
tems with lower Nhin will have statistically lower Asngi 
values and higher Aaii values, while the opposite holds for 
systems with Nbin > 5. This behavior is simply related 
to Poissonian statistics. In general, the derived Asngi and 
Aaii values corroborate the picture revealed by Figure 4, 
i.e., the upper limits on the observed variations in individ- 
ual DLA are between 0.2 and 0.4 dex. In several systems 
- especially those with higher S/N - the upper limits are 
below 0.2 dex at the 95% c.l. and Abest < 0.1 dex. Within 
these systems, it is difficult to accommodate any signifi- 
cant variation in the intrinsic abundances. 

4. IMPLICATIONS 

4.1. Nucleosynthetic Enrichment and Dust Depletion 

The results presented in the previous section demon- 
strate that the majority of damped Lya systems have rel- 
ative ionic column abundances which are uniform to better 
than 0.2 dex. In fact, the majority of components show 
departures < 0.1 dex. This analysis reveals an important 
characteristic of high z galaxies: the gas 'clouds' which 
comprise these protogalaxies apparently have very similar 
physical properties and nucleosynthetic enrichment histo- 
ries. 

In several respects this uniformity is unexpected. Con- 
sider first the variations one might predict from nucleosyn- 
thetic enrichment. We stress that the metallicities of the 
individual gas 'clouds' comprising a damped Lya system 



are unknown because the HI line-profiles are unresolved. 
Variations in metallicity greater than a factor of 10 are pos- 
sible if not probable. The components may have very differ- 
ent enrichment levels and could be expected to show differ- 
ent nucleosynthetic patterns. The magnitude of variations 
in the nucleosynthetic patterns, however, would depend 
on the actual range of metallicities and the detailed star 
formation histories. While extremely metal-poor Galac- 
tic halo stars exhibit significant variations in their rel- 
ative abundances (e.g. Co/Ni, Zn/Fc: McWilliam et al. 
1995; Johnson 2002), the majority of Galactic stars with 
[Fe/H] = —2 to —1 dex have nearly uniform relative abun- 
dances (Gratton & Sneden 1991; Fulbright 2000). There 
are, however, a number of stars and globular clusters in 
this metallicity range which show significantly lower a/Fe 
ratios (e.g. Ru 106, Pal 12, lo Cen, etc.; Brown, Waller- 
stein, & Zucker 1997; Nissen & Schuster 1997; Smith et al. 
2000; Stephens & Boesgaard 2002). Although these stars 
may have a distinct origin from the bulk of the Galactic 
halo stars (e.g. accreted from other galaxies), they reveal 
nucleosynthetic variations at all metallicity. 

This point is further emphasized by the differences be- 
tween the Milky Way abundances and other local galaxies 
(e.g. LMC, SMC, d'Sph galaxies: Venn 1998; Shetrone, 
Cote, & Sargent 2001). In the SMC, stars with metal- 
licity [Fe/H] = -0.7 have solar Si/Fe ratios (Venn 1998) 
in stark contrast with values of [Si/Fe] > -1-0.2 observed 
at the same metallicity in the Milky Way. Similarly, the 
few measurements of dwarf spheroidal galaxies (Shetrone, 
Cote, & Sargent 2001; Bonifacio et al. 2000; Shetrone 
2002) show a dispersion in Si/Fe of 0.25 dex at metallicity 
between 1/10 and 1/100 solar. Altogether, one observes a 
range in a/Fe of at least 0.3 dex both within a given galaxy 
and from galaxy to galaxy. One might have predicted a 
similar dispersion in the majority of high z protogalaxies. 
On the other hand, it is important to note that Kobulnicky 
& Skillman (1997) found only small deviations in the HII 
regions throughout NGC 1569. They concluded that the 
by-products of star formation within this galaxy reside in 
a hot 10^ K phase which is unmixed with the general ISM. 
If this process is important at high z, it would impact the 
abundance variations within DLA. 



Table 2 
DUST SUMMARY 



Gal 


Phase" 


D(Zn) 


D(Si)'^ 


D(Fc) 


[Zn/Fc] 


[Si/Fc] 


Ref 


MW 


warm halo 


0.0 


-0..3 


-0.6 


+0.6 


+0..3 


1 


MW 


warm disk 


0.0 


0.4 


1.4 


+1.4 


+1.0 


1 


MW 


cold disk 


-0.1 


-1.3 


-2.2 


+2.1 


+0.9 


1 


LMC 


warm halo 


0.0 


-0.2 


-0.4 


+0.4 


+0.2 


3 


LMC 


warm disk 


0.0 


-0.1 


-0.9 


+0.9 


+0.8 


3 


LMC 


cold disk 


0.0 


-0.4 


-1.2 


+1.2 


+0.8 


3 


SMC 


warm halo 


0.0 


0.0 


-0.6 


+0.6 


+0.6 


2 


SMC 


warm disk 


0.0 


0.0 


-1.4 


+1.4 


+1.4 


2 


SMC 


cold disk 


0.0 


0.0 


-2.2 


+2.2 


+2.2 


2 



"Crude phase description, largely based on depletion level 
''Depletion factor (e.g. [X/Zn]) 

References. — Key to References - 1: Savage and Sembach (1996); 
2: Welty et al. (1999); 3: Welty et al. (2001) 
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Fig. 5. — Same as Figure 3 except the deviations are plotted as a function of DLA metallicity. For comparison in the lower panel, we show 
the abundance variations measured in the LMC (at [M/H] = —0.3) and SMC (at [M/H] = —1) for Si/Fe (squares) and Zn/Fe (stars). 



In addition to the implications for nucleosynthetic en- 
richment in the damped systems, our observations impact 
the nature and prevalence of dust in high z protogalax- 
ies. Table 2 summarizes the typical depletion levels and 
gas-phase abundance ratios for various interstellar phases 
of the IMilky Way, LJVIC, and SMC (Savage and Sem- 
bach 1996; Welty et al. 1997, 1999, 2001). Because the 
[Si/Fe] and [Zn/Fe] ratios vary significantly from phase 
to phase, most sightlines penetrating these galaxies show 
variations in excess of 0.5 dex. For Si/Fe, the situation 
is complicated because Si is refractory, e.g., one observes 
roughly the same Si/Fe ratio in the JVIilky Way warm and 
cold disk gas even though the depletion levels differ by 
> 0.5 dex (Table 2). Nevertheless, most sightlines through 
the Galaxy and JVlagellanic Clouds show large variations 
in Si/Fe. This is in stark contrast with our observations 
of the damped Lya systems. 

When comparing the damped systems against these galax- 
ies, it is important to keep in mind that even the SMC is 
more metal-rich than most damped systems. If dust for- 
mation is very sensitive to gas metallicity, one might ex- 



pect variations in differential depletion for only the most 
metal-rich damped systems. Interestingly, we observe no 
correlation between the magnitude of abundance varia- 
tions and metallicity. Figure 5 plots the deviations from 
uniformity and their significance as a function of the gas 
metallicity (either Si/H, Zn/H or S/H) including a set of 
values for the LMC and SMC (Welty et al. 1999, 2001). 
Note that the systems with the highest metallicity have 
among the most uniform abundances in our sample. It 
appears that the ISM of the majority of DLA - indepen- 
dent of metallicity - is significantly different from mod- 
ern galaxies. This begs an obvious but puzzling question: 
How is a galaxy like the SMC different from the majority 
of DLA such that it regularly shows highly depleted gas? 
We also emphasize that the low values of Asngi indicate it 
is rare in the damped systems to have even a single, highly 
depleted 'cloud'. If highly depleted gas exists at high red- 
shift, its filling factor is very small in most protogalaxies. 

Finally, independent of expectations from abundance 
variations observed in local galaxies, there is a long stand- 
ing debate on whether the damped Lya abundances reveal 
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an underlying Type la or Type II SN enrichment pattern 
(e.g. Lu ct al. 1996; Vladilo 1998; Prochaska & Wolfe 
2002; Vladilo 2002). These two patterns have a/Fe ratios 
differing by « 0.3 dex and imply very different enrichment 
histories (e.g. Lindner, Fritzc-von Alvcnslebcn, & Frickc 
1999). The source of this debate is the difficulty in dis- 
entangling the nucleosynthetic pattern from the depletion 
patterns of the gas-phase abundances. Although this mat- 
ter remains unsettled, it is possible - if not likely - that 
some damped systems arc dominated by Type la SN en- 
richment with solar a/Fe ratios while others exhibit the 
a-enriched patterns associated with Type II SN. Most rel- 
evant to the current discussion, however, is the prospect 
that both processes contribute to the abundances of indi- 
vidual protogalaxies. The observed chemical uniformity of 
the damped Lya systems places tight constraints on pos- 
sible nucleosynthetic variations. For example, scenarios 
where the majority of gas has nearly solar Si/Fe or Zn/Fe 
ratios with a single component enhanced by 0.3 dex are 
ruled out in over half of the DLA. 

4.2. The Physical Nature of the DLA 

To fully appreciate the implications of the DLA abun- 
dance uniformity, one would like to interpret the results in 
terms of the physical nature of the damped systems. Un- 
fortunately, the morphology of the damped Lya systems 
remains an open and difficult question. While the low- 
ion kinematics are remarkably well described by a thick, 
rotating disk (Prochaska & Wolfe 1997b, 1998), the ob- 
servations also permit scenarios within the CDM hierar- 
chical cosmology that describe DLA as multiple merging 
'clumps' bound to individual dark matter halos (Haehnelt 
et al. 1998; McDonald & Miralda-Escude 1999; Mailer et 
al. 2001). Furthermore, several authors have also claimed 
that the damped Lya systems might be explained by out- 
flows from SN winds (Nulsen, Barcons, & Fabian 1998; 
Schaye 2001). The difficulty in distinguishing between 
these scenarios stems from the limited information pro- 
vided by the 'pinhole' nature of quasar absorption line 
studies. Although a statistical sample of systems can con- 
fidently rule out many proposed scenarios (Prochaska & 
Wolfe 1997b), a given sightline cannot uniquely express 
the galaxy's morphology and even a sample of sightlines is 
limited. In the following, we reflect on several of the lead- 
ing scenarios but avoid outflow scenarios because these 
models have only been crudely developed and it remains 
unclear whether they can succcssfiilly account for the ve- 
locity fields of most DLA. Nevertheless, many of the con- 
clusions for the other models apply to the wind scenarios. 

Assume, first, that the DLA primarily arise from gas 
'clouds' within a single, gravitationally bound system. In 
terms of nucleosynthetic variations, the observed unifor- 
mity might be easily accommodated within this scenario. 
In particular, one can assume that the majority of gas and 
stars comprising the protogalaxy have had a similar chem- 
ical enrichment history. There might be regions of local 
star formation which would exhibit a dispersion in nucle- 
osynthetic enrichment, but these regions would have small 
filling factor and the by-products may not be directly dis- 
persed into the general ISM (e.g. Kobulnicky & Skillman 
1997). For larger, more massive systems one might require 
that the system is well mixed, i.e., the dynamical time 



is short compared to the age of the system. At z ~ 2, 
the galactic dynamical time is t^yn < 50 Myr assuming 
tdyn ^ '^Rd/vc and Rd = Ar2oo/\/2 where A is the spin pa- 
rameter and r2oo is the virial radius (Mo, Mao, & White 
1998). Therefore, it is reasonable to expect that the nucle- 
osynthetic enrichment within a single protogalaxy would 
not lead to variations greater than 0.1 dex and possibly 
much less. 

In contrast with nucleosynthetic variations, the absence 
of differential depletion variations is striking. The low 
^sngi values for most damped systems indicate it is rare 
for sightlines to penetrate even a single, highly depleted 
gas 'cloud'. This is in stark contrast to the LMC, SMC, 
and Milky Way where one routinely intersects clouds with 
a broad range of depletion levels. It appears that the 
damped systems are comprised almost exclusively of ei- 
ther undepleted gas or gas with warm halo depletion pat- 
terns. If gas with high levels of depletion exist at high z. its 
porosity is very small, at least within the surface density 
contour defining a damped Lya system. Perhaps highly 
depleted gas exists only in the inner, densest regions of a 
protogalaxy which has very small cross-section relative to 
the outer layers of gas. This observational constraint may 
have important implications for the detailed nature of star 
formation in young galaxies. 

Within the CDM cosmogony, the scenarios which best 
reproduce the damped Lya kinematics of larger At; sys- 
tems are those where quasar sightlines penetrate multi- 
ple 'clumps' or satellites within a given dark matter halo. 
These satellites represent the remnants of mergers between 
the present dark matter halo and a previous (now con- 
sumed) halo. In this scenario, the uniformity of abundance 
ratios within the damped Lya abundances constrains the 
relative abundances of all of the satellites comprising a 
damped system. Presumably, these satellites have unique 
enrichment histories and unique physical characteristics, 
e.g., gas fraction, molecular gas content, density, velocity 
fields. In terms of nucleosynthesis, this implies the satel- 
lites share a similar enrichment history. 

On theoretical grounds, the observed a/Fc ratio corre- 
sponds to the ratio of Type II to Type la SN and, therefore, 
crudely assesses the rate of star formation (e.g. Tinsley 
1979). As described above, galaxies in the local universe 
exhibit a range of a /Fe ratios at the range of metallicities 
observed for the damped Lya systems. Altogether, one ob- 
serves a variation in a/Fe of »0.3 dex both within a given 
galaxy and from galaxy to galaxy. The accreted satellites 
which comprise a damped system in these CDM scenarios 
might be expected to have a similar range of metallicity, 
unique star formation histories, and therefore exhibit a 
range in a/Fe. The constraints our observations place on 
these nucleosynthetic histories depend on the distribution 
of a/Fe values. The observations are inconsistent at the 
95% C.I., for example, with a bimodal distribution of a/Fe 
ratios with values of and -1-0.3 dex. On the other hand, 
only a few profiles are inconsistent with a scenario where 
all of satellites have a/Fe values ranging uniformly from 
to +0.3 dex (i.e. Ac,;/ < 0.15). We stress, however, that the 
majority of Af,est values are < 0.1 dex indicating that even 
this distribution is not favored.. Higher S/N data might 
lower the 95% c.l. limits to 0.1 dex and pose a true chal- 
lenge to multiple clump scenarios, especially as chemical 
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evolution models aehieve greater sophistication. 

Similar to the 'single system' models, the constraints on 
dust depletion variations are strong. With the exception 
of a few damped systems (notably Q133H-17, Q0347-38), 
the observations require that the 'clumps' have very simi- 
lar differential depletion. In other words, most of the gas 
comprising the individual satellites in these CDM models 
must have a very low porosity of gas with highly depleted 
levels. From the results in previous papers (e.g. Pettini 
ct al. 1997; Prochaska & Wolfe 2002), wc knew the in- 
tegrated depletion levels in the damped systems is both 
small and uniform. We now appreciate that the com- 
ponents comprising individual damped systems also are 
largely undepleted and with very similar differential de- 
pletion. 

5. SUMMARY AND CONCLUDING REMARKS 

We have performed a quantitative analysis of chemi- 
cal abundance variations along the sightlines through 13 
damped Lya systems. These systems each have veloc- 
ity width At; > 40 km s^-'^, i.e., large enough base-line to 
examine chemical uniformity. We performed several sets 
of Monte-Carlo simulations to place upper limits on rela- 
tive abundance variations. To our surprise, the majority 
of DLA exhibit a high degree of chemical uniformity. In 
most cases, the dispersion in intrinsic abundanc;es the 
combined effects of all abundance variations - is less than 
0.2 dex and the best values are generally < 0.1 dex. 

For the ratios examined in this paper, there arc three 
potential sources of intrinsic abundance variation: (1) nu- 
cleosynthetic enrichment; (2) differential depletion; and 
(3) photoionization corrections. In terms of depletion, the 
observed DLA uniformity contrasts with variations in dif- 
ferential depletion observed along sightlines through the 
SMC, LMC, and Milky Way. Although the latter galaxies 
have higher metallicity than most DLA, the difference in 
depletion variations is not explained solely by differences 
in metallicity; there is no correlation between abundance 
variations and metallicity in the DLA (Figure 5). For an 
unknown reason or reasons (perhaps related to dust de- 
struction mechanisms), the nature of dust at high z is 
qualitatively different from the local universe. Considering 
nucleosynthetic processes, the DLA abundances are more 
uniform than the dispersion in nucleosynthetic enrichment 
of the Milky Way as traced by stellar abundances. Fur- 
thermore, one observes a greater dispersion between stars 
in various galaxies within the local group than that ob- 
served in the gas of the DLA. Finally, photoionization cal- 
culations suggest that Si+ZFe"*" can vary by over 0.2 dex 
between a neutral gas and a highly ionized gas (Howk 
& Sembach 1999; Vladilo et al. 2001; Prochaska et al. 
2002a) . Our observations suggest that the cross-section of 
HII regions or similarly ionized gas is small for galaxies 
with Af(HI) > 2 X 10^0 cm-2. 

The chemical uniformity of the DLA poses an impor- 
tant constraint on the nature of high z protogalaxies. We 
have argued that 'single-system' scenarios might reason- 
ably account for the observed abundance invariance. Cen- 
tral to this conclusion, however, is that these protogalax- 
ies have a small filling factor of highly depleted gas. If 
star formation is linked to dust-depleted molecular clouds, 
then our results indicate the gas relevant to star forma- 



tion encompasses a very small cross-section. This conclu- 
sion is consistent with the low fraction of DLA showing 
molecular gas (Petitjean et al. 2000). In contrast with 
the single-system models, we contend the observed uni- 
formity presents a unique challenge to the multiple-clump 
scenarios favored by CDM. These protogalactic clumps or 
'satellites' do not share a common gas reservoir and should 
have unique physical characteristics (density, metallicity) 
and enrichment histories. It remains to be demonstrated 
whether these clumps might express vc;ry similar differen- 
tial depletion and nucleosynthetic enrichment patterns. 

Before concluding, we wish to comment on several fu- 
ture observational efforts which could improve upon this 
paper. First, we emphasize that the results for many of the 
DLA in the current sample are limited by signal-to-noise. 
Follow-up observations of the brighter quasars might re- 
veal uniformity at the 0.01 dex level and place qualitatively 
tighter constraints on processes of metal production and 
dust formation. Similarly, higher resolution data would 
enable an investigation in velocity bins of a few kms~^. 
Second, we would like to repeat our analysis using Zn 
and Si separately to better isolate the effects of deple- 
tion and nucleosynthesis. Variations in the Zn/Si ratio 
have been observed in at least one low z DLA (Pettini et 
al. 2000) and we would like to examine similar variations 
at high z. Third, it is also important to examine abun- 
dance ratios which relate to physical processes separate 
from the Type la vs. Type II enrichment and differential 
depletion emphasized in this paper. For example, a study 
on the variations of the N/a ratio within DLA would im- 
pact our understanding of star formation timescales and 
possibly the universality of the initial mass function (e.g. 
Prochaska et al. 2002b). Similarly, because the Ar'^/Si"'" 
or 0°/Si+ ratios are particularly sensitive to photoioniza- 
tion, an analysis of these ions would help reveal processes 
of ionization within the DLA. Finally, a comparison of the 
C II* fine-structure line with Si and Fe resonance transi- 
tions bears on the nature of the protogalactic ISM and 
ultimately star formation rates (Wolfe, Prochaska, & Ga- 
wiser 2002). Together these observations would help re- 
veal the detailed physical properties of the gas comprising 
DLA protogalaxies. 
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